Scaling relation of the anomalous Hall effect in (Ga,Mn)As 
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We present magnetotransport studies performed on an extended set of (Ga,Mn)As samples at 4.2 
K with longitudinal conductivities a xx ranging from the low- to the high-conductivity regime. The 

1 is extracted from the measured longitudinal and Hall resistivities, 
i — l 



anomalous Hall conductivity a xy 
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regime to a scaling relation a x 
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a xx for low-conductivity samples is observed. This scaling 
relation is consistent with a recently developed unified theory of the anomalous Hall effect in the 
framework of the Keldysh formalism. It turns out to be independent of crystallographic orientation, 
growth conditions, Mn concentration, and strain, and can therefore be considered universal for 
low-conductivity (Ga,Mn)As. The relation plays a crucial role when deriving values of the hole 
concentration from magnetotransport measurements in low-conductivity (Ga,Mn)As. In addition, 
the hole diffusion constants for the high-conductivity samples are determined from the measured 
longitudinal conductivities. 

PACS numbers: 75.50.Pp, 72.20.My 
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I. INTRODUCTION 

In III-V-based diluted magnetic semiconductors such 
as (Ga,Mn)As, the ferromagnetic interaction between the 
localized Mn spins is mediated by mobile holes^ Spin 
polarization and spin-orbit coupling of the carriers re- 
sult in a large anomalous Hall effect (AHE), which has 
been discussed intensively and controversially since the 
19<S0s£d£JLLIL The total Hall resistivity of ferromagnets 
originates from the sum of the normal Hall effect and the 
AHE 9 



Pyx 



p (NH) 
ryx 



Rah VoM, 



(1) 



where Rq and Rah are the ordinary and the anomalous 
Hall coefficients, respectively. B denotes the magnetic 
induction and M the magnetization perpendicular to the 
layer. Until a few years ago, the AHE was exclusively 
ascribed to scattering anisotropies induced by spin-orbit 
interaction? 9 , known as the extrinsic skew scattering 6 and 
side jumpS mechanisms. These mechanisms lead to con- 
tributions to Pyx 11 "* showing a linear and quadratic de- 
pendence on the longitudinal resistivity p xx , respectively. 
Accordingly, Eq. |T]) can be rewritten as 



Pyx = RqB + C\p,r + fj/', 



(2) 



where c\ and ci are constants for magnetic fields high 
enough to saturate M or, to be more precise, to sat- 
urate the hole spin polarization. An anomalous Hall 
term, quadratic in p xx , is also obtained for the in- 
trinsic (scattering independent) mechanism related to 
the Berry phased It was first associated with the 
AHE in (Ga,Mn)As by Jungwirth et al. in their sem- 
inal paper Ref. Illl where they supposed the anoma- 
lous Hall conductivity in metallic (Ga,Mn)As to orig- 
inate from strongly spin orbit-coupled holes acquiring 
a Berry phase during their motion on the spin-split 



Fermi surfaces. Meanwhile the Berry phase mechanism 
is considered the dominant contribution to the AHE 
in metallic (GajMnJAs^ii 2 -' 1 ^' 1 ^ Whereas most pa- 
pers concerning the AHE in (Ga,Mn)As are primarily 
focused on the metallic regim e 11 ' 12 ! 13 ] 14 ! 15 ] 16 , only a few 
systematic studies concentrate on the low-conductivity 
regimei 13 ! 17 ' 18 Allen et alJ£ examined the AHE in dig- 
itally doped (Ga,Mn)As structures with hopping trans- 
port on the insulating side of the metal insulator transi- 
tion (MIT). Recently, Shen et alJ^ focused on the low- 
conductivity regime of (Ga,Mn)As with p xx between 0.01 
and 2 il cm. They experimentally ascertained a scaling 
relation p^ H ' ) cx p1 x with a scaling exponent of 7 = 0.5. 
The same value for 7 was obtained in Ref. [l^ for a high- 
resistivity (Ga,Mn)As sample showing large magnetore- 
sistance. Neither in Ref. [l6| nor in Ref. [lT] a theoreti- 
cal explanation for this value was presented. In Ref. [H, 
the authors analyzed a series of (Ga,Mn)As samples and 
observed a crossover of the anomalous Hall coefficient 
Rah cx p xx from the low-resistivity regime dominated by 
the Berry phase with n — 2 to the high-resistivity region 
with n = 1.3. 

In this work, we experimentally find a scaling expo- 
nent of 7 = 0.4 in the low-conductivity regime, which 
is consistent with a unified theory of the AHE, based 
on the Keldysh formalism, recently developed for multi- 
band ferromagnetic metals by Onoda et a L 19 ' 20 This 
model predicts three regimes depending on the car- 
rier scattering time: In the "clean" limit (extremely 
high-conductivity samples), the extrinsic skew scatter- 
ing mechanism with Py^ H ^ cx p xx dominates. This 
regime is beyond the conductivities achievable in the 
current (Ga,Mn)As growth technology. In the interme- 
diate regime (high- conductivity samples), the intrinsic 
Berry phase mechanism with p\^^ cx p 2 xx prevails, and 
in the "dirty" limit (low-conductivity samples) the intrin- 
sic contribution is strongly damped, resulting in a scaling 
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relation p yx oc p x '%. As shown below, the conductivities 
in the (Ga,Mn)As samples under study correspond to the 
intermediate and dirty regime. Therefore, according to 
the model of Onoda et al., Eq. ([SJ has to be replaced by 

p yx = R B + c 2 p 2 xx + c z p x i . (3) 

Whereas in nonmagnetic p-type semiconductors the hole 
density p can be easily derived from the ordinary Hall co- 
efficient Rq — 1/ ep, obtained by magnetotransport mea- 
surements at low magnetic fields, the determination of p 
in (Ga,Mn)As is complicated by the dominant anomalous 

contribution p^^ in Eq. {!]). To overcome this problem, 
p yx and p xx are usually recorded as a function of B above 
the saturation field and Rq is derived by fitting Eq. ([1]) 
to the experimental values of p yx . In general, its value 
is significantly affected by the choice of the relation be- 
tween />yz H ' ) and p xx . In Refs. and [20L the AHE is 
described in terms of the anomalous Hall conductivity 
ciy 1 ^ and longitudinal conductivity a xx instead of Py^^ 
and p xx . Since in (Ga,Mn)As p yx <C p xx has been exper- 
imentally established, the conductivities and resistivities 
are related by 



III. RESULTS AND DISCUSSION 

In Fig. [TJ the Curie temperature Tc, estimated from 
the maximum of the temperature dependent longitudi- 
nal resistivity p xx r^ and the magnetoresistance (MR) 
Apxx{B)/p xx (0T) at 5=14.5 T of the samples under 
study are plotted as a function of the zero-field longi- 
tudinal conductivity a xx [B = T). While Tc increases 
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II. EXPERIMENTAL DETAILS 

In order to investigate the scaling relation between 
ciy H ^ and a xx in (Ga,Mn)As, we have studied the AHE 
in a set of twenty samples with longitudinal conduc- 
tivities covering about one order of magnitude from 
39 cm -1 to 442 fi -1 cm -1 . The samples with thick- 
nesses from 30 to 250 nm and Mn concentrations between 
2% and 5% were grown by low-temperature molecular- 
beam epitaxy. Samples of different quality were ob- 
tained by variation of the substrate temperature and/or 
the V/III flux ratio, as well as by using undoped semi- 
insulating (001)- and (113) A GaAs wafers or partially 
relaxed (001) (Ga,In)As/GaAs templates. All layers 
showed the typical (Ga,Mn)As reflection high energy 
electron diffraction (RHEED) pattern during the growth 
with no visible second phase of hexagonal MnAs clus- 
ters. Further experimental details concerning the sample 
growth can be found in Refs. 22,23,24. The samples with 
high conductivities a xx >180 cm -1 were obtained by 
postgrowth annealing for 1 h at 250 °C in air. Through- 
out this paper, samples with a xx > 150 fi^cm -1 at 
4.2 K are referred to as metallic. 1000 /imx 300 pm Hall 
bar structures were prepared from the samples by stan- 
dard photolithography and wet chemical etching. The 
transverse and longitudinal resistances were measured by 
magnetotransport at 4.2 K applying magnetic fields up to 
14.5 T. 



FIG. 1: Curie temperature Tc (open circles) and absolute 
value of MR at 14.5 T (closed circles) of the samples under 
study. The vertical dashed line represents a rough boundary 
between the high- and low-conductivity regime. 

with increasing a xx , the absolute value of the MR de- 
creases and drops below 7% in the metallic samples. 
Figures and [3] show the measured Hall resistivity p yx 
and the MR of two representative samples A and B, 
covering the cases of low (a xx w 57 cm -1 ) and 
high longitudinal conductivity (a xx w 442 cm -1 ) at 
14.5 T, respectively. The increase of p yx and p xx below 
~1 T stems from the anisotropic magnetoresistance. 22 
The Hall resistivity of sample A shown in Fig. [DJa) is 
dominated by the large negative MR, which can be at- 
tributed to the suppression of weak localization effects in 
the presence of a magnetic field B. 25 In order to evalu- 
ate the applicability of Eqs. |J5J| and both equations 
were used to fit the experimental Hall data between 2 
and 14.5 T. The inset of Fig. Ufa) clearly shows that 
the measured curve is better reproduced using Eq. © 
(solid line), with C3p xx being the dominant term, than 
using Eq. © (dashed line). Furthermore, the dashed 
curve underestimates the measurement below 2.5 T, un- 
physically implying that M increases with decreasing 
B. The fit using Eq. ([3]) yields a hole concentration of 
p = 1.8 x 10 20 cm~ 3 , in contrast to the less reasonable 
low value of p — 7 x 10 19 cm~ 3 obtained from Eq. 
Note that the same values are obtained if the fit is started 
at magnetic fields higher than 2 T. Comparison between 
Fig. [2(b) and GJb) shows that the MR for B > 2 T in 
the metallic sample B is about two orders of magnitude 
smaller than in sample A, resulting in an extremely weak 
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FIG. 2: (Color online) (a) Measured Hall resistivity of the 
low-conductivity sample A (open circles) and fit curves using 
Eq. ([2} (dashed line) and Eq. ([3} (solid line). The inset shows 
the Hall resistivity and the fits on an enlarged scale, (b) 
Magnetoresistance of sample A. 

dependence of p xx on B. Therefore, the AHE makes only 
a marginal contribution to the slope of p yx and, as a con- 
sequence, Eqs. @ and ^ yield the same fit curves as 
demonstrated in the inset of Fig. [3Ja). In both cases, 
a hole concentration of p = 8.5 x 10 20 cm -3 is derived. 
Generally, we find that in metallic samples with a suf- 
ficiently small MR, fits to p yx yield the same values for 

Pyx^ and p irrespective of the exponents used for p xx 
making it impossible to separate the different contribu- 
tions to the AHE. Therefore, Eq. ^ can be used to model 
the p yx curves for the whole set of (Ga,Mn)As samples 
under study. 

Figure Q] illustrates the main finding of this paper. It 

shows the anomalous Hall conductivity ci^ H ' ) as a func- 
tion of g xx at 14.5 T, determined from the measured 
resistivities using Eqs. and dH). Irrespective of 
the growth conditions, Mn concentration, and substrate 
orientation, the data for the low-conductivity samples 
(&xx < 100 il -1 cm -1 ) follow the same universal scaling 
relation 

with s = 0.0076 O - 6 cm - 6 and a = 2 - 7 = 1.6. 
This value of a is in agreement with the theoretical pre- 
dictions made in Refs. [l^ and for multiband ferro- 
magnetic materials in the "dirty" limit. As mentioned 
in Sec. U the data presented by Shen et al*il for low- 
conductivity (Ga,Mn)As follows a very similar scaling 

relation a x ^ H ^ oc cr^ 5 (dashed-dotted line), which they 
claim to be incompatible with existing scattering theo- 
ries. If the MR of a sample is large, the scaling relation 



FIG. 3: (Color online) (a) Measured Hall resistivity of the 
high-conductivity sample B (open circles) and fit curves using 
Eq. ([2]) (dashed line) and Eq. <[3j) (solid line). In this case the 
fits according to Eqs. ((2| and ([3| fall together. The inset 
shows the Hall resistivity and the fits on an enlarged scale, 
(b) Magnetoresistance of sample B. 
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FIG. 4: (Color online) Anomalous Hall conductivity ai y as 
a function of the longitudinal conductivity a xx at 14.5 T and 
4.2 K. The solid line represents a universal scaling <riy H ' oc 
ai f. The dashed dotted line indicates the scaling given in Ref. 
IT7I . The plateau value of a^y = 20 cm -1 is attributed 
to the Berry phase effect. 



in Eq. ([5]) can also be tested by plotting a xy '(B) as a 
function of the corresponding magnetic-field dependent 
longitudinal conductivity a xx (B). Figure [5] shows the 
appropriate log-log plot for the low-conductivity sample 
A which exhibits the largest MR(14.5 T)=-42% of the 
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samples under study. We obtain a = 1.63 (7 = 0.37) tudinal charge transport. In Fig. EJ^a), a xx measured at 

B = T is plotted as a function of p. From the data 
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FIG. 5: (Color online) log- log plot of the anomalous Hall con- 
ductivity cr±^ of the low-conductivity sample A as a function 
of the magnetoconductivity a xx (B) at 4.2 K. The solid line 

represents a fit yielding ai^j^ = 0.0076 a xx z . 

in good agreement with the aforementioned theory and 
the value 7 = 0.5 derived by Pu et alJ& from a low- 
conductivity (Ga,Mn)As sample. Interestingly, even for 
the insulating (Ga,Mn)As samples with hopping trans- 
port as studied in Ref. |l_8| the scaling exponent 7 turned 
out to be close to 0.4. Due to the fact that the scal- 
ing exponent a — 1.6, or equivalently 7 = 0.4, has 
also been found in numerous other low-conductivity fer- 
romagnets, independent of the details of the underlying 
transport process (see e .g. Refs. [3 ll9ll20l . and references 
therein, and Refs. we are led to the conclu- 

sion that it has to be considered universal and applies 
to low-conductivity (Ga,Mn)As irrespective of the spe- 
cific material parameters. For samples with longitudi- 
nal conductivities between ~100 and ^200 f2 _1 cm , a 
gradually decreases to a — with increasing a xx in qual- 
itative agreement with the theory presented in Refs. [19| 
and I20L For the high-conductivity samples with a xx > 
200 cm -1 , a xy becomes independent oia xx amount- 
ing to o'i^ H ' ) = 20 f2 _1 cm -1 . This value is in quantita- 
tive agreement with theoretical data in Ref. [TU calculated 
within the framework of the Berry phase theory of the 
AHE for a disordered sample with p = 6.0 x 10 20 cm~ 3 . 
Since the mean value of p derived for the metallic samples 
under study is also close to 6.0 x 10 20 cm -3 (see below), 
our result for the high-conductivity regime corroborates 
the intrinsic origin of the AHE in metallic (Ga,Mn)As. 
As demonstrated in Fig. 2J a significant linear contribu- 
tion (riy H ' ) oc a xx due to skew scattering can be ruled out 
for the investigated range of a xx . This is in agreement 
with the calculations of Nagaosa et al. 4 giving an esti- 
mate a xx ;§> 1000 f2 _1 cm _1 for the onset of dominant 
skew scattering mechanism. 

In the remainder of this work we focus on the longi- 



FIG. 6: (Color online) (a) Longitudinal conductivity a xx at 
T as a function of the hole concentration p. The regres- 
sion line indicates the proportionality DhNh(Ep) oc p in the 
metallic regime, (b) Open circles represent the diffusion con- 
stant of the holes derived from Eq.([6j|. The value extracted 
by Neumaier et al.— is indicated by an open square. 

we expect the MIT to occur between p ss 1 x 10 20 cm~ 3 
and p ps 2 x 10 20 cm -3 . The simple relation a xx = ep/ih 
yields a hole mobility of \ih — 3 cm 2 V _1 s — x , inferred 
from the slope of the regression line (straight line) for 
the metallic samples, in satisfactory agreement with the 
values given in Refs. [3(3 and [3l|. In fact, the adequate 
description of the hole transport in metallic (Ga,Mn)As 
at 4.2 K, neclecting small quantum corrections, is given 
by the generalized Einstein relatio n 32 ! 33 

a xx = e 2 D h N h (E F ) , (6) 

where Dh denotes the hole diffusion constant and 
Nh(Ep) the density of states (DOS) at the Fermi en- 
ergy Ep. Thus, the experimentally observed linear in- 
crease of a xx with increasing p in the metallic regime 
indicates a proportionality DhNh{Ep) oc p. Only few 
paper a 34 ' 35 have been published so far concerning the 
hole diffusion constants in (Ga,Mn)As in one- and two- 
dimensional systems, and only little is known about the 
values of Dh in bulk (Ga,Mn)As. Using Eq. ([6]) and cal- 
culating Nh(Ep) in the framework of a 6 x 6 k ■ p Zener 
modelf 2 - we have determined the hole diffusion constant 
Dh of the metallic samples as a function of p. Figure [Sfb) 
shows that Dh varies between ~ 3.0 x 10~ 5 m 2 s _1 and 
~ 6.5 x 10 -5 m 2 s _1 for the metallic samples. The ex- 
trapolation of these values to lower hole concentrations 
agrees with the one found by Neumaier et ali^ The val- 
ues Nh{Ep) used in Eq.© are shown in Fig. E[a). The 
absolute value of the Fermi energy, extracted from the 
same calculations, is displayed in Fig. [T^b). The ex- 
change splitting parameter Bq , primarily influencing the 
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covering about one order of magnitude. For the samples 
in the low-conductivity regime, a universal scaling 



p(10 20 cm" 3 ) 



FIG. 7: (a) DOS N h (E F ) at the Fermi energy in (Ga,Mn)As 
calculated by means of the Zener model of Dietl et al.— The 
value experimentally derived from Neumaier et al>2£ is indi- 
cated by an open square, (b) Absolute value of the Fermi 
energy Ef calculated as a function of the hole density p. 



or equivalently, pl^^ 



0.4 

XX 



relation a x ^ H ^ oc o xx 
found, which can be modeled within a fully quantum- 
mechanical transport theory for multiband systems 
recently developed in Refs. [H and In the metallic 
regime, the scaling changes to a scattering independent 
anomalous Hall conductivity criy H ' ) = 20 il" 1 cm" 1 

(Pyx^ = 20 f2 _1 cm -1 x p 2 IX ) due to the Berry phase 
effect. The skew scattering mechanism can be ruled out 
to play a significant role for the range of a xx under study. 
Consequently, Eq. ^ instead of Eq. (T2]) , commonly cited 
in the literature, has to be used for the determination 
of the hole concentration in high- and low-conductivity 
(Ga,Mn)As from high-field magnetotransport data. 
In addition, hole diffusion constants between 
3.0 x 10" 5 m 2 s -1 and 6.5 x 10" 5 m 2 s" 1 for the 
metallic samples are obtained using the measured values 
of a xx and calculating the corresponding density of 
states Nh(Ep) at the Fermi energy within the Zener 
model. 



values of Nh(Ep) and Ep, was taken proportional to 
p as suggested by the study in Ref. 23. The calcula- 
tions yield an increase of \Ep | and Nh(Ep) from ^180 



- 2.6 xlO 46 J' 
respectively. 



L m 3 with increasing p > 3.5 x 10 2 



IV. SUMMARY 
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